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We study the effect of QCD corrections to tfie gg ^>- H ^>- W^W~ signal-background interference 
at the LHC for a heavy Higgs boson. We construct a soft-coUinear approximation to the NLO and 
NNLO corrections for the background process, which is exactly known only at LO. We estimate its 
accuracy by constructing and comparing the same approximation to the exact result for the signal 
process, which is known up to NNLO, and we conclude that we can describe the signal-background 
interference to better than 0(10%) accuracy. We show that our result implies that, in practice, a 
fairly good approximation to higher-order QCD corrections to the interference may also be obtained 
by rescaling the known LO result by a /^-factor computed using the signal process. 



I. INTRODUCTION 

Search for the Higgs boson at the LHC has been a 
remarkable success so far. Indeed, both the ATLAS 
and CMS collaborations have announced the discovery 
of a new boson, whose properties are compatible with 
that of the Standard Model Higgs particle, with mass 
mh ~ 125 GeV. Both collaborations also excluded ad- 
ditional Higgs-like bosons in a large mass range rrih ^ 
600 GeV [l|, |2|. The interpretation of the excesses ob- 
served in various production and decay channels, as orig- 
inating from a single spin-zero particle, was made possi- 
ble by detailed theoretical predictions for the Higgs boson 
production and decay rates, sec Ref. [3| for an overview. 

However, these experimental results do not imply that 
there are no additional Higgs-like bosons with masses 
600 GeV < mh < 1 TeV. In fact, the search for such 
particles is well underway [J] . In the Standard Model, as 
the Higgs boson becomes heavier, its total decay width 



grows rapidly Th 



li thanks to contributions of the 



longitudinal electroweak bosons: for m/j ~ 600 GeV, the 
width is close to 120 GeV. Since the finite-width effects 
change the distribution of the invariant masses of the de- 
cay products of the Higgs boson, their understanding is 
important for developing experimental search strategies. 

There are two finite width effects that influence the 
Higgs boson lineshape. First, the Higgs propagator must 
assume the Breit-Wigner form in the resonant regime 
l/(s — m^) — )• l/(s — m^ -I- imfiTji). While this modifica- 
tion is literally correct for a light (and therefore narrow) 
Higgs boson, for a heavy Higgs, it must be modified; the 
proper way to do this was subject to a significant discus- 
sion in recent literature, see Refs. |5|, |6| and references 
therein. The second effect is the interference with the 
background. Note that, in principle, the two effects are 
not completely independent of each other since modifica- 
tions of the Breit-Wigner form for the propagator change 
the very definition of the "background" in the resonance 



region, but discussion of these subtleties is beyond the 
scope of this paper. 

Our goal is to consider the interference of the signal 
process (75 — ^ _ff — > W^W~ and the background process 
gg — > W~^W~ for a heavy Higgs boson^. This interfer- 
ence was first computed at leading order in Refs. 0, Q. 
Although the gg — >■ W^W~ amplitude appears at one 
loop, it is enhanced at the LHC by the large gluon flux, 
making the interference effects non- negligible. An obvi- 
ous shortcoming of Refs. 0,13 is that their analysis of the 
interference is performed at leading order in perturbative 
QCD as far as the Higgs boson signal is concerned. This 
is unfortunate since, for the Higgs boson signal, higher 
order QCD corrections are extremely important, as they 
enhance the total rate by more than a factor two |9l4l]| . 
It is therefore interesting to explore their impact on the 
signal-background interference. 

Such an endeavor, however, is highly non-trivial. In- 
deed, a full NLO and NNLO QCD calculation of back- 
ground amplitudes requires evaluation of two- and three- 
loop 2 — )■ 2 Feynman diagrams which is beyond the reach 
of the current computational technology. On the other 
hand, it is well-known [l2| that for the Higgs boson sig- 
nal a large fraction of radiative corrections is captured 
by the soft-collinear approximation. Since this approxi- 
mation should be particularly suitable for the description 
of a heavy Higgs boson, we construct a soft-collinear ap- 
proximation for the entire gg — > W~^W~ amplitude that 
includes both the signal and the background and study 
the impact of these corrections on the interference. 

This paper is organized as follows. In Section [H] we 
sketch the construction of the soft-collinear approxima- 
tion. In Section IIIII we present numerical results. We 
conclude in Section HVl 



^ For the light m^ = 125 GeV Higgs boson the interference is 
negligible if proper signal-selection criteria are applied Ql . 



II. SETUP 

We begin by describing the setup of our computation. 
We are interested in higher order QCD corrections to 
the interference between the signal process gg — >■ -ff — )■ 
W~^W~ and the pure QCD background gg — !• W^W' . 
We compute these corrections in the soft gluon approx- 
imation, which is known to describe the full NLO and 
NNLO Higgs cross section to very good accuracy. We 
will numerically assess the accuracy of our approxima- 
tion in Sec. IIIII by comparing it with known NLO and 
NNLO results for the signal process. 

The cross section for the production of a W^W~ pair 
with invariant mass Q^, fully differential in the kinemat- 
ics variables of the two M^'s, is given by 

dcr(T,2/, {fli},^^) = / da;ida:2dz/g(.Ti,^F)/g(a;2,AiF) 
X b{T -x-YX2z)Aa[ z,y,{6'j},as,^-,^- ) (1) 

V Mr Mf / 

where jg is the gluon distribution, and da is the differ- 
ential partonic cross section for the process 

5(pi) +5(P2) -> W^[vw^) -t- W-[vw-)+ X, (2) 

with {pw+ +Pw-Y ~ Q^\ Mf and /Ir are the factoriza- 
tion and the renormalization scales, ots = asinn) is the 
strong coupling constant at the scale /Ir, t = Q^/s. We 
denote by y the rapidity of the W pair, and by {9i} a 
generic set of variables describing the kinematics of the 
decay products of the W~^W~ system in the hadronic 
center-of-mass frame; they are related to the correspond- 
ing variables y, {9i} in the partonic center-of-mass frames 
by a boost with rapidity ycm = 5 1^ |J- , and thus the 9i 
are functions of {0i}, xi, X2 and z. 

In the soft (z — >■ 1) limit, the rapidity distribution 
of the W^W~ pair is entirely determined by the inclu- 
sive cross section [l3l - [l5| . up to corrections suppressed 
by powers of (1 — z), and the partonic cross section in 
Eq. ^ takes the form 



dff 



Mr Mf 



da(")({^,;},a.)zGfz,a„^,^ 



Mr M 



(3) 



where da^^\{9i} , as)S{l~ z) is the leading order partonic 
cross section, and G[z, as,Q^/n^, Q^/Mf) is the inclusive 
coefficient function computed in the soft limit, i.e. (up 
to the explicit z factor) the inclusive partonic cross sec- 
tion normalized to the leading order in such a way that 
G(z,«,)=<5(l-z) + 0(a,). 

In the same limit, the momenta of the W bosons in 
the partonic centcr-of-mass frame are given by 



with 9 the W boson scattering angle in the partonic 
center-of-mass frame, and /3 = y^l — Arn^Y/Q"^ (for sim- 
plicity, wc have assumed that the W^-bosons are on-shell, 
but we will not make this assumption in the sequel) . The 
kinematics of the process in the soft limit is therefore the 
same as the leading order kinematics, except that the 
total energy squared is rescaled by a factor z. 

The boost that relates the partonic and hadronic 
center-of-mass frames is fixed by taking for the momenta 
of the colliding gluons either pi = zxiPi, p2 = X2P2 
or pi ~ XiPi, P2 = ZX2P2, where P12 are four- 
momenta of the colliding protons [1J|- Alternatively, 
one may also take as momenta of the colliding gluons 
Pi = ^/zxiPi, p2 = \fzx2P2 [131 • These two choices 
coincide in the soft limit up to terms suppressed by two 
powers of (\ — z) [l5| and, in fact, give very similar results 
for observables considered in this paper. We will make 
the first choice at NLO, where it is actually exact, while 
at NNLO we will take the average of the results obtained 
with either choice cases. 

We now turn to the explicit form of the coefficient func- 
tion G'(z,as,Q^/MR,Q^/MF)j which contains the core of 
our soft-coUinear approximation. We first sketch the im- 
portant features of the soft gluon approximation and its 
modifications by focusing on the next-to-leading order. 
Further details on this, including re quir ed modifications 
at NNLO, can be found in Refs. [li,li3- 

Working to NLO accuracy and in the soft limit and 
neglecting all non-singular terms, we write the function G 
as (we suppress explicit scale dependence for simplicity) 



G(z,a,) = 5(l-z) 



27r 



27r2 



%CaDx (z) + ( — Ca + ci ) 5(1 - z) 



(5) 



where Viiz) = [ln'(l — z)/(l — z)] and ci is the ratio 
of the infrared regulated higher-order virtual contribu- 
tions to the cross section and the leading order cross 
section for gg — > W^W^ , see [ia| for its proper defi- 
nition.^ For our purposes, the important feature of this 
formula is that non universal NLO corrections for the 
process gg — > WW only enter through the coefficient Ci . 
This is because only emissions from external gluon lines 
in each diagram contribute to the amplitude in the soft 
limit. For the signal-only proces s gg —>■ H —>■ WW, c\ 
is known both in the infinite mt [ISl . Il9| approximation 
and for finite rrit |20| . The determination of c\ for the 
interference would require the evaluation of complicated 
gg — ?► W~^W~ amplitudes which is beyond existing tech- 
nical capabilities. 

However, we note that the value of ci can be ob- 
tained without any computation in the kinematic limit 



Pw± 



I,±(3sm9,0,±l3cos9 



(4) 



^ Because we consider here the 2 — > 2 scattering process, ci does 
depend on the scattering angle. We assume that this depen- 
dence is mild and systematically ignore it in this paper. Partial 
justification for this assumption is given below. 



4m'^ 



In this limit, the in- 



< Q^ < 4mj, rUb ^ rut. 
terference is dominated by the contribution of longitudi- 
nally polarized W bosons, which can be obtained from 
QCD corrections to the production of two neutral scalars 
gg — > HH in the heavy top mass limit |21j. Since both 
the box contribution for gg — > HH and the triangle con- 
tribution for gg ^ H are described by the same effec- 
tive Lagrangian, the virtual QCD corrections should be 
identical in the two cases. Although the assumptions 
Q^ ^ 4?7i( , mt ~ rrib are not really justified, we take the 
value for ci that is obtained in that limit as a reference 
value, and estimate the sensitivity of the final result to 
its variations. 

The soft approximation of Eq. ^ is of course only de- 
fined up to sublcading terms. An optimal choice of sub- 
leading terms can be found 17[ by using a combination 
of analiticity arguments in Mellin space, and information 
on universal subleading terms in the z — >■ 1 limit, aris- 



ing partly from the exact soft-gluon kinematics 15 1 and 
partly from universal collinear splitting kernels 



12 



A discussion of this optimal soft approximation is be- 



yond the scope of this paper, and we refer to Ref. [17 1 
for a full discussion. Here, we note that the best approx- 
imation proposed in [l3| (called soft2 there) effectively 
amounts to performing in Eq. ([5]) the replacement 



V,{z)~^V,{z) + 5V,{z), 



dV,{z) = {2-3z + 2z 



In'i^ 

2\ Vz 



1 



In'(l-z) 
1-z 



(6) 



where 6'Di{z) is an ordinary function (not a distribution). 
In what follows, we will call the approximation based on 
Eq. ([5]) with such replacement a "soft-collinear" approxi- 
mation. We will quantify the impact of subleading effects 
by comparing this improved soft-collinear approximation 
to a purely soft result. 

At higher orders the soft approximation Eq. ^ is also 
known: sec e.g. Eq. (79) in [id. We improve it anal- 
ogously to Eq. (l6|), see Ref. [13] for details. This soft- 
collinear approximation is the basis for the NLO and 
NNLO numerical results for the signal and the interfer- 
ence that we discuss in the next Section. 



III. NUMERICAL RESULTS 

We consider the process gg -^ W^{e^i^)W^{e^D) at 
the LHC for two values of the centcr-of-mass energy: 
yi = 8 TeV and ^/s = 13 TeV. We take the Higgs mass 
to be mh = 600 GeV, and its total decay width to be 
Th ~ 122.5 GeV [2j]. All numerical results presented be- 
low are obtained with a fixed-width Breit-Wigner func- 
tion. We have checked that use of the running-width in 
the Breit-Wigner propagator [2J| leads to results for the 
signal and interferences that differ by an amount that 
is below our accuracy goal, and we expect that same is 
likely to be the case for a full treatment of finite-width 
effects [SJ, jg. Moreover, we have found that the QCD 



TeV 



13 TeV 



NLO 



NNLO 



NLO 



NNLO 



exact 

soft-collinear 

iV-soft 



2.150 
2.187 
2.135 



2.78 

2.820 

2.700 



2.074 2.67 
2.127 2.730 
2.073 2.607 



TABLE I: ii'-f actors for the inclusive Higgs-only cross section 
in the narrow width approximation, with m^ = 600 GeV, 
computed using the exact theory, our best soft-collinear ap- 
proximation, and an unimproved soft approximation (see text 
for details). The (N)NLO result is computed using (N)NLO 
PDFs, while the reference LO cross section is always com- 
puted with NLO PDFs. Numerical results are obtained using 
the code H. 



radiative corrections are insensitive to the propagator, to 
the accuracy we work to. We let both the M^-bosons de- 
cay Icptonically and reconstruct all kinematic variables 
from the charged lepton and neutrino momenta. We take 
the W total width to be F^y ~ 2.11 GeV and heavy quark 
masses mt = 172.5 GeV and m\, ~ 4.4 GeV. 

We use the NNPDF2.3 PDF set at NLO and 
NNLO, with as{niz) — 0.118. Throughout this paper, 
we set the rcnormalization and factorization scales equal 
to the Higgs boson mass /j,r = ^p = "i^. In constructing 
our soft-collinear approximation, we retain the exact mt 
and mh dependence where available. For example, we 
use the exact value of ci, Eq. ([5]), for the signal process, 
while for the analogous 0{al) coefficient C2 we use the 
value computed in the infinite m^t (point-like) approxi- 
mation. Note that with this choice, all logarithmic terms 
at NNLO have the exact mt and mi, dependence, while 
the coefficient of the (5(1 — z) term is only approximate. 
As mentioned in Sect. [Hi for the interference we take the 



result in the 



'w 



< Q^ < mf, mb 



mt 



limit as our 



reference value. 

To assess the quality of the soft-collinear approxima- 
tion, we first test it against the signal-only gg ^ H pro- 
cess at NLO and NNLO. Results are shown in Tab. HI for 
two values of the collider energy. The iC-factors com- 
puted (without including the Higgs decay) using the ex- 
act theory'^ are compared to those obtained with our soft- 
collinear approximation, or with the so-called A'^-soft ap- 
proximation, defined in Ref. [17[. The latter amounts to 



approximating the partonic cross section with the inverse 
Mellin transform of a pure A'^-space soft approximation. 



^ At NNLO, an exact result valid for large Higgs masses is not 
currently available. For our result, we use the exact result at 
NLO [lS| plus the point-like result at O(a^), improving it with 
those mt, mi, dependent terms which are fully determined by 
lower orders (which include all soft-collinear terms). We have 
checked that the result obtained in this way is stable upon varia- 
tion of small- 2: terms up to the accuracy shown in Table IT] which 
is a consequence of the dominance of soft-collinear terms for a 
heavy Higgs boson at the LHC [2711 . 



STeV 



13 TeV 



8TeV 



Vs = 13 TeV 



LO NLO NNLO LO NLO 



NNLO 



GH 0.909 1.99(5) 2.6(1) 
aH^ 1.188 2.6(1) 3.4(3) 

LO 



„ /„Ll 

„ /„LO 
(THi/Oui 



2.19(5) 2.8(1) 
2.2(1) 2.9(2) 



3.77 8.1(2) 10.3(5) 

4.56 9.7(4) 12.5(9) 

— 2.14(5) 2.7(1) 

— 2.13(9) 2.8(2) 



TABLE II: Results (in fb) for the Higgs-only cross sec- 
tion GH and the signal+interference cross section a Hi, with 
mh = 600 GeV. No cuts on the final state applied. The errors 
represent the uncertainty on the soft-coUinear approximation 
and on the unknown background coefficients, estimated as 
explained in the text. 



in which only powers of In N and constant terms are kept. 

Both approximations reproduce the exact result to 
0(3%) or better in all configurations. At yG = 8 TeV, 
where the soft-coUinear terms are expected to domi- 
nate [27| , our soft-coUinear approximation reproduces the 
exact result to better than 0(2%), while at higher energy, 
y/s = 13 TeV, the agreement deteriorates slightly, be- 
cause non-soft terms become relatively more important. 
However, whereas at NNLO the soft-coUinear approxi- 
mation is more accurate than the A^-soft, at NLO the 
opposite happens. This occurs because numerically the 
A'^-soft approximation happens to be closer to the exact 
result than our improved soft-coUinear one in the small- 
N limit. Since the small- A^ limit is beyond the region of 
applicability for both of these approximations, we con- 
sider this feature to be accidental but note that one can 
improve both of these approximations by matching them 
to the correct small- iV limit [231 • In what follows we use 
the soft-coUinear approximation as the default and take 
the spread of values between the soft-coUinear and the 
A^-soft approximations as an estimate of the uncertainty 
due to deficiencies of these approximations in the small- 
N region. 

We have also checked the reliability of our approxi- 
mation for differential distributions when decays are in- 
cluded. Indeed, at NLO accuracy, we find that our ap- 
proximate results for the lepton pt and rapidity distribu- 
tions and for the lepton invariant mass mu distribution 
are in good agreement with the full result obtained from 
MCFM [H]. 

Having assessed the accuracy of our approximation, we 
can now apply it to study higher order corrections to the 
signal-background interference. As explained in the pre- 
vious Section, we need the exact leading order prediction 
for the interference. We extract it from Ref. [7|, as im- 
plemented in MCFM. For the Higgs boson signal, we use 
the exact expression obtained as discussed above. For the 
background, we include the contributions of all the three 
quark generations, see [3] for details. We also need the 
infrared- regulated virtual cross section ci , and the analo- 
gous NNLO coefficient C2. As already mentioned, we take 
the signal values for these coefficients ci.2 as a reference, 





LO 


NLO 


NNLO 


GH 


0.379 


0.83(2) 


1.07(5) 


GHi 


0.427 


0.93(3) 


1.20(7) 


'^h/g^P 


— 


2.19(5) 


2.8(1) 


f^ Hi /GHi 


— 


2.19(7) 


2.8(2) 



LO NLO 



NNLO 



1.55 3.29(8) 4.2(2) 

1.66 3.5(1) 4.5(2) 

— 2.13(5) 2.7(1) 

— 2.12(6) 2.7(1) 



TABLE III: Same as Table HIl but with Higgs-based cuts on 
the final state. See text for details. 



and study the impact of virtual corrections on the inter- 
ference by varying ci_2 in the range — 5ci^2 < ci^2 < 5ci_2- 

We first discuss the impact of QCD corrections on 
the inclusive cross section. Following Ref. Q, we 
compare the signal-only cross section uh with the 
background-subtracted cross section aui = o'gg^ww ~ 
agg-fww\hg only, which includes interference effects. We 
report our results for the signal only cross section an and 
the signal+interference cross section ant for ci^2 = ci^2 
in Tablcim To facilitate the comparison with the results 
of Ref. [7|, LO results are computed using NLO PDFs. 
For the signal, the quoted error is obtained by comparing 
our soft-coUinear approximation to the A^-soft approxi- 
mation. For the background, we also consider the ad- 
ditional uncertainty coming from independently varying 
the ci_2 coefficients for the first two and the third gener- 
ation in the — 5ci^2 < ci_2 < 5ci^2 range. This leads to an 
uncertainty of about 6% on the interference predictions 
which, combined with the uncertainty of the soft approx- 
imation, gives an overall uncertainty of about 8 — 9% at 
NNLO, see Table [Hi This uncertainty is of same order 
of magnitude as the current uncertainties in the Higgs 
production rate ctnnlo related to higher-order QCD ra- 
diative corrections, PDF and a^ uncertainties etc, see 0]. 
We conclude that our approach to estimate higher order 
corrections to the signal-background interference in the 
Higgs production offers a robust framework and adequate 
phenomenological precision. 

Wc turn to a discussion of the impact of the interfer- 
ence in a more realistic setup, by imposing selection cuts 
on leptons and neutrinos. Apart from the standard ac- 
ceptance cuts on the lepton rapidity r]i , lepton transverse 
momentum pt and missing energy ^^ , 

|r;,|<2.5, pt>25GeV, ^^ > 20 GcV (7) 

we impose additional signal-enhancement cuts, linearly 
extrapolating numerical values given in Ref. |3C| . To this 
end, we require at least one lepton with pt > 130 GeV, 
and impose the following cuts on the lepton invariant 
mass mil , aziniuthal separation A(j)ii of the two leptons 
and transverse mass of the W'^W^ pair m^: 

mil < 500 GeV, A0/, < 3.05, 

120 GeV <m^<mh. (8) 

Wc note that we have validated the soft-coUinear ap- 
proximation at NLO QCD against MCFM for the differ- 
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FIG. 1: Lepton azimuthal distance Ai^;; distribution in the fully inclusive case (left pane) and with experimental cuts (right 
pane) computed with the NNLO QCD soft-collinear approximation described in the text. Dots show the rescaled MCFM result 



for the signal dcr^SD ^ 



X i^NNLo/^NLO, where ii'(N)NLO is the inclusive iC"- factor. 
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FIG. 2: Same as Fig. [U but for the lepton invariant mass mu distribution. 



ential distributions, so that we believe that our results 
are reliable even when cuts on the final state are im- 
posed. We report our results in Tab. IIIII We see that 
the impact of the interference is mildly (but notably) re- 
duced when the Higgs-selection cuts are applied to the 
final state particles. Note also that radiative corrections 
to the interference are rather similar to corrections to the 
signal cross section. 

We conclude this Section by showing the effect of the 
interference on selected kinematic distributions at the 
13 TcV LHC. In Fig. [1] we plot the difference of the 
azimuthal angle Ac/);/ of the two charged leptons with 
(right pane) and without (left pane) Higgs-selection cuts. 
In Fig. [5] we do the same for the invariant mass of the 
charged leptons mu- We plot the NNLO QCD results ob- 
tained with our soft-collinear approximation as described 
in Scct.ini using ci^2 = ci^2 for the interference case. We 
see that the Higgs-sclcction cuts reduce the importance 
of the interference, as already seen in the total rate. 

An interesting feature of our results is that our approx- 
imation reproduces, to a good accuracy, all the kinematic 
distributions as obtained with MCFM. In particular, all 



the distributions can be perfectly reproduced by rescal- 
ing the MCFM leading order distributions by the inclu- 
sive NNLO i^-factor. For the signal, we also compare 
our NNLO approximation against the known NLO distri- 
butions, rescaled by the NNLO/NLO inclusive iiT-factor 
(also shown in the plots). Also in this case, the agree- 
ment is excellent; the only exception is the azimuthal 
angle distribution where differences are seen at large rel- 
ative angles. This is due to the fact that our soft-collinear 
approximation does not reproduce the effects of a hard 
emission, which modify the angular distribution. Note, 
however, that the azimuthal angle cut plays an insignif- 
icant role in separating the heavy Higgs boson from the 
background so that the impact of this mismatch on cor- 
rections to the interference is minor. 



IV. CONCLUSIONS 

We have estimated the impact of QCD radiative cor- 
rections on the signal-background interference in gg — > 
H —J- W'^W~ process for a heavy Higgs boson. We con- 



structed a soft-collincar approximation to higher-order 
QCD corrections and verified its vahdity by comparing it 
to exact results for gg ^f H , including kinematic distri- 
butions of the Higgs decay products. We find that QCD 
radiative corrections enhance the signal-background in- 
terference by a significant amount which, however, is very 
similar to the perturbative QCD enhancement of the sig- 
nal cross section. 
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